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The electronic and optical properties of self-assembled InN/GaN quantum dots (QDs) are inves- 
tigated by means of a tight-binding model combined with configuration interaction calculations. 
Tight-binding single particle wave functions are used as a basis for computing Coulomb and dipole 
matrix elements. Within this framework, we analyze multi-exciton emission spectra for two different 
sizes of a lens-shaped InN/GaN QD with wurtzite crystal structure. The impact of the symmetry 
of the involved electron and hole one-particle states on the optical spectra is discussed in detail. 
Furthermore we show how the characteristic features of the spectra can be interpreted using a sim- 
plified Hamiltonian which provides analytical results for the interacting multi-exciton complexes. 
We predict a vanishing exciton and biexciton ground state emission for small lens-shaped InN / GaN 
QDs. For larger systems we report a bright ground state emission but with drastically reduced 
oscillator strengths caused by the quantum confined Stark effect. 

PACS numbers: 78.67.Hc, 73.22.Dj, 71.35.-y 
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I. INTRODUCTION 

The great topical interest in semiconductor nanostruc- 
tures is not only based on the realization of some of the 
paradigm of elementary quantum mechanics^, but also 
by the wide range of possible applications, ranging from 
new- 3 -^ and extremely efficient light sources^ to building 
blocks for quantum information technology 7 -*^. Further 
motivation is provided by the possibility to study the 
effect of reduced dimensions on carrier transport^ and 
optical propertie d 1 ^ 12 ' 13 . 

The confinement of the carriers in all three dimensions 
on a nanometer scale is achieved, for example, by em- 
bedding regions of a material with a smaller bandgap in 
a matrix of a wider bandgap material. A widespread 
method of creating such quantum dots (QDs) is the 
Stranski-Krastanow growth mod o 14 : 15 ' 16 . 

In this fast evolving research field, nanostructures 
based on conventional group-Ill nitrides like A1N, GaN, 
and InN gained more interest in recent years. Compared 
to group-Ill arsenide semiconductor materials, nitride- 
based nanostructures have the advantage that it is possi- 
ble to span a much larger spectral range, presently from 
amber to ultraviolet, by properly alloying together the 
three building blocks and thereby engineering the direct 
band gap of these materials^ 7 .. 

Nevertheless and despite the intense research on ni- 
tride QDs over the last decade the understanding of this 
material system is - compared to other III-V materials - 
still in its infancy. 

From the theoretical point of view, one challenge is 
the proper inclusion of effects that stem from the altered 
atomistic structure of the underlying lattice. While most 
of the nitrides can crystallize both in the zinc-blende and 
the wurtzite phase, the latter is by far more stable^. 
Additionally, the strong built-in field needs to be consid- 
ered for, the mostly applied growth along the c-axis. In 



contrast to many other III-V semiconductors, the spin- 
orbit coupling in the nitrides is rather weaki 7 - so that the 
calculation of the electronic states in terms of a simple 
effective-mass approximation is not possible due to strong 
valence band mixing effects. Instead, especially for small 
nanostructures, a microscopic description of the single- 
particle states based on, for example, a tight-binding 
model or pseudo-potential calculation is necessary. 

In the past many experimental and theoretical inves- 
tigations addressed the optical properties of III-V and 
II- VI QD structures based on, e. g., InGaAs/GaAs or 
CdSe/ZnSe. One important aspect concerned the analy- 
sis of absorption and emission spectra of QDs as a func- 
tion of the excitation densityi^^^^^ 2 ^ 3 * 2 ^. A 
common result is that the optical processes mainly in- 
volve 'diagonal' transitions, that are connecting, for ex- 
ample, s-shell electrons with s-shell holes or p-shell elec- 
trons with p-shell holes. In envelope-function approxi- 
mation this is traced back to dipole matrix elements cal- 
culated from the Bloch functions, and interband tran- 
sition amplitudes determined by the product of dipole 
matrix elements and the overlap of the electron and hole 
envelope functions. This picture has been proven to be 
very fruitful for conventional III-V materials and is often- 
times carried over to the nitride system where emission 
spectra are then calculated using these 'diagonal' exci- 
tonsj 26 ' 27 : 28 : 29 : 30 : 31 Our analysis shows, however, that for 
nitride QDs deviations from this picture are possible. We 
identify situations, in which the emission is dominated 
by recombination of 'skew' excitons, that are excitons 
which consist of s-shell electrons and p-shell holes or vice 
versa. This led to the prediction of dark exciton and 
biexction ground state emission. 32 Such a modified se- 
lection rule cannot be explained in an simple effective- 
mass picture but requires the inclusion of band-mixing 
effects in a way that properly accounts for the symmetry 
of the underlying atomic lattice. In addition to the im- 



2 



portant changes of the optical properties involving one 
and two excitons, the transitions that occur in systems 
with a higher population of excitons differ strongly from 
those known from the InGaAs/GaAs system. While for 
smaller nitride QDs the "skew" excitons are predicted, 
we find that in large nitride QDs the energetic order of 
the s-shell and p-shell for holes is interchanged, with the 
p-shell being lowest in energy. This leads to optically 
active exciton transitions. 33 However, as a result of the 
quantum confined Stark effect (QCSE), with increasing 
height of the QDs the corresponding dipole matrix ele- 
ments decrease drasticall y 28 i 33 ' 34 due to the separation of 
the electron and hole wave functions in the strong inter- 
nal electric field. In combination with the dark exciton 
state in small nitride QDs, this makes the application of 
nitride QDs more challenging. 

In this work we study the multi-exciton emission spec- 
tra in nitride QDs and discuss the resulting complicated 
peak structure in detail. Restricting ourselves to the two 
lowest shells allows even a semi-analytic description of 
the problem. The single-particle states are deduced from 
an atomistic tight-binding model as briefly outlined in 
Sec.|n] Even without a detailed calculation of the single- 
particle states, the symmetry considerations of Sec. IIIII 
can be used to draw conclusions regarding the dipolc- 
matrix elements and the multi-exciton spectra. In Sec. lIVI 
the Hamiltonian of the interacting charge carrier system 
and the configuration interaction (CI) scheme are pre- 
sented. Results for an initial filling with up to six excitons 
are provided in Section IV Bl and major trends are dis- 
cussed in Section [V Cl in terms of an approximate Hamil- 
tonian. Further details of the multi-exciton spectra are 
analyzed in Section IVT1 In the subsequent Section IVlII the 
influence of the strong built-in field is investigated. Fi- 
nally, the spectra for a larger nitride QD are presented in 
Section lVIIII and compared to those found for the smaller 
structure. 



II. ELECTRONIC STRUCTURE 

For a proper treatment of the single particle states 
in nitride systems, an atomistic multiband approach, 
like pseudopotential or tight-binding (TB) models, is re- 
quired. General aspects for the developement of such a 
TB approach and the calculation of dipole and Coulomb 
matrix elements are discussed in detail in Ref. 33. In the 
following Section III Ai we briefly summarize the main 
ingredients of the subsequently used TB model, while 
Sections III Bl and III CI are devoted to details of the cal- 
culation of single-particle states and interaction matrix 
elements. 



A. Tight-Binding model and built-in field 

For the calculation of the single-particle states we use 
a TB model with an sp 3 basis |q;)r, that contains one 
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FIG. 1: (color online). Top view of the QD structrue with 
the first two bound shells for electrons (left) and holes (right). 
Depicted are isosurfaces of the probability density with 20% 
(blue) and 80% (red) of the maximum value are shown. Only 
one state is visualized for the two-fold degenerate p-shell, as 
the other looks alike. The corresponding energies (E^' 2 3 ) are 
measured relative to the valence band maximum of the bulk 
GaN. 



s-state (a = s) and three p-states (a — p x ,p y ,p z ) per 
spin direction at each atomic site R. The i th TB single- 
particle wave function is given by 

|&) =J>>r<4(R) 

with the expansion coefficients denoted by c^(R). We 
assume a lens-shaped InN QD, grown in (OOOl)-direction 
on top of an InN wetting layer (WL), that is embedded 
in a GaN matrix. For the numerical evaluation, a finite 
wurtzite lattice within a supercell with fixed boundary 
conditions is chosen. The small spin-orbit coupling and 
crystal-field splitting are neglected^. Furthermore, we 
neglect the strain induced displacement of the atoms. For 
the used QD geometry, a more realistic inclusion of strain 
effects does not change the symmetry of the system and 
henceforth the general statements discussed in this paper 
should not be affected. In the following we consider two 
different QD sizes with diameters d = 4.5, 5.7 nm and 
heights h — 1.6,2.3 nm, respectively. In both cases a 
WL thickness of one lattice constant c is assumed. 

In contrast to cubic III-V semiconductor heterostruc- 
tures, e.g., InAs or GaAs, the III-V wurtzite nitrides 
exhibit considerable larger built-in fields 35 , which can 
significantly modify both the electronic and the optical 
properties. This field enters via the electrostatic poten- 
tial as a site-diagonal contribution V^(R) = — e0 p (R) to 
the TB Hamiltonian, 3 ^! 3 ! 



B. Single-particle states 

The small QD (d = 4.5 nm, h = 1.6 nm) confines three 
bound states for the electrons. The corresponding prob- 
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ability densities of the first two bound shells for electrons 
and holes are presented in Fig. [T] In the calculation, the 
influence of the built-in field was included. The dominant 
contribution for the electron single-particle states stems 
from the atomic s-orbitals, while for the hole states a 
strong intermixing of the atomic p-orbitals is observed. 
This is indicative for the fact that a suitable description 
of the one-particle states and the resulting optical prop- 
erties requires a multi-band approach and cannot be ac- 
counted for in a one-band effective mass approximation. 

If one compares the single-particle states with and 
without the inclusion of the built-in field, 33 one finds 
that, in the presence of the field, the electron states are 
squeezed into the cap of the QD while the hole states are 
constraint to a few atomic layers at the bottom near the 
WL. This effect is even more pronounced for the larger 
QD (d = 5.7 nm, h = 2.3 nm). In this case a clear 
spatial separation of electron and hole wave functions is 
observed, which lowers the direct spatial overlap between 
the two and leads to strongly reduced dipole matrix ele- 
ments. Furthermore, as known from the QCSE, the built- 
in field leads to an overall red-shift in the single-particle 
transition energies. In addition, we find for the large QD, 
that the ordering of the energy levels is affected by the 
built-in field. Specifically we obtain that in the presence 
of the built-in field the level ordering for the holes is re- 
versed so that the ground state is formed by the two-fold 
degenerate states 4>2 3 (p-shell in Fig. [I]) while the first ex- 
cited state 4>i (s-shell in Fig. [TJ) is nondegenerate. This is 
in contrast to the observed level ordering of the large QD 
in the absence of the built-in field or of the small QD with 
and without the field. In these cases, the ground state is 
nondegenerate (4>\) and the first two excited states are 
degenerate (^2,3) ■ Such a behavior has been reported be- 
fore for other Q D sy stems and a detailed discussion can 
be found in Ref. [33l 



Evaluation of dipole- and interaction matrix 
elements 



Based on the TB single-particle wave-functions one 
can determine dipole- and Coulomb matrix elements that 
are crucial for the calculation of optical properties. As 
emphasized above, a TB model represents an atomistic 
approach to describe the electronic structure of low- 
dimensional heterostructures. Explicit knowledge about 
a basis set of localized (Wannier) states is not required 
for the calculation of one-particle energies and wave func- 
tions, because the Hamiltonian matrix elements between 
the different orbitals are treated as parameters within the 
TB model. What enters the TB calculation are only the 
basic assumptions about the localized (atomic) orbitals: 
symmetry, spatial orientation 33 , and orthogonality. 

For the calculation of interaction matrix elements one 
needs - at least in principle - the localized basis states. 
For the Coulomb matrix elements, however, the explicit 
knowledge of the atomic orbitals is not required in prac- 



tice, as these matrix elements are dominated by long- 
range contributions in which the local orbitals act as 
point charges. The structure of the localized orbitals is of 
significance only for on-site and nearest-neighbor interac- 
tions, which in our calculation contribute less than 5 % to 
the total Coulomb matrix elements. These findings are in 
agreement with Ref. 39. Thus, the matrix elements are 
approximated by a sum over the TB coefficients c^(R) 
at atom sites R, R' with orbital indices a, (3: 



V ljM =Jd 3 rJ d 3 r' ^(r)^(r')^(r-r')0 fc (r')^(r) 

« E E c(R)4*( R ')4( R >L(R) ^( R - R ') 

RR' ccP 

e 2 

with VCR - R ) = : r for R ^ R 

V. I A _ IT) T}/ I 



4 7 r£ £r|R-R / 



and V(0) = 4- / 



d 3 rd 3 r'- -. 

47T£ £r|r - r'| 



The labels i,j,k,l refer either to electron or to hole 
states in order to consider the repulsive electron-electron 
and hole-hole interaction, as well as the attractive 
electron-hole interaction. The considerably smaller ma- 
trix elements of the electron-hole exchange interaction 
are neglected. The calculation of the onsite integrals 
V(0) involves the integration over the volume of the unit 
cell v uc and can be done quasi-analytically by expansion 
of the Coulomb interaction in terms of spherical harmon- 
ics^ The electronic charge and the vacuum dielectric 
constant are denoted by e and eo, respectively. We use 
the InN dielectric constant e r = 8.4 according to Ref. |4l| 
since the wave functions are almost completely confined 
inside the QD. 

For the calculation of dipole matrix elements d^j 1 cx 
{<j)f |r|0j), the explicit structure of the localized orbitals 
is needed as the dipole-operator is dominated by short- 
range contributions. Though being not orthogonal at 
different sites, standard Slater orbitals 4 ^, have been used 
in earlier calculation a 43 ' 44 within orthogonal TB models. 
While they include the correct symmetry properties, the 
missing orthogonality limits their applicability. To over- 
come this problem, we use numerically orthogonalized 
Slater orbitals 3 - 3 - and account for the slight non-locality 
of the dipole operator— and the underlying anion-cation 
structure of the crystal by including also contributions 
from up to second-nearest neighbors. 



III. SYMMETRY CONSIDERATIONS 

The specific symmetry of the system under considera- 
tion plays an important role in the prediction of energy 
degeneracies and optical selection rules. The overall sym- 
metry of the problem is determined by two factors: (i) 
the crystal symmetry of the underlying lattice and (ii) 
the symmetry of the QD geometry or, more generally, 
the geometry of the heterostructure4!i 
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The interplay between the dot geometry and the un- 
derlying lattice is well illustrated by the example of the 
lens-shaped QD. The geometry of the QD has a Coc V 
symmetry^ In the effective mass approximation or in 
the k • p theory the underlying atomic structure is not 
resolved, so that one obtains degenerate energies for the 
lens-shaped QD i 46 ' 48 ' 49 However, if one takes the crys- 
tal structure into account, as it is done, in a T B 33 i 50 
or pseudo-potential calculation 5 - 1 -, the symmetry will be 
reduced and degeneracies can be lifted.— For a lens- 
shaped QD grown on a wurtzite lattice along the (0001)- 
direction, the symmetry is reduced to C$ v . If a QD with 
the same geometry is grown in the zinc-blende phase in 
(00 Indirection, one is left with an even lower Ci v symme- 
try. As we will discuss in the next subsection, the former 
group is still rich enough to predict degeneracies, while 
the latter is too small to support degenerate eigenstates. 



A. Energy spectrum 

The C^v group is generated by the rotation around the 
z-axis with an angle of 27r/3, denoted by C^ Zl together 
with a mirror reflection at a plane perpendicular to the y- 
axis, denoted by ICi y ^ These operations commute with 
the Hamiltonian, but not among themselves. From this 
non-abelian character of the group one can immediately 
conclude that there must exist energetically degenerate 
states. In the case of the wurtzite QD under considera- 
tion, one example of such degenerate states are the p-shell 
levels given in Fig. [T] 

This degeneracy in the wurtzite structure is especially 
worth mentioning in view of a recent paper— in which it 
is argued that a description at the atomic level in zinc- 
blende-based QDs should remove the p-shell degeneracy. 
The argument is essentially the following: the four-fold 
rotation is part of the zinc-blende point group only if 
followed by the inversion, but many high symmetry (lens- 
shaped, pyramidal) QDs forbid the inversion and with it 
the Ca v symmetry. One is therefore left with a Czv point 
group, which is abelian and therefore insufficient for p- 
degeneracy. 

This conclusion is correct, but it cannot be carried over 
to the wurtzite QDs, because in wurtzite crystals the 
three-fold rotation holds without inversion. If the QD 
geometry is sufficiently symmetric (lens-shaped, hexago- 
nal pyramids) then the problem retains the C 3v symme- 
try, which leads to p-degeneracy. In papers predicting p- 
shell splitting for QDs in wurtzite lattices and neglecting 
the spin-orbit coupling 53 ! 54 , either boundary conditions 
or discretization meshes spoil the correct C 3v symmetry 
of the problem. Therefore the small p-shell splitting ob- 
tained there is a numerical artefact. On the contrary 
the splitting induced by the inclusion of the spin-orbit 
coupling 5 -, albeit quite small, is physically correct. 

To summarize the discussion, for a lens-shaped QD and 
in the absence of the spin-orbit coupling a degeneracy of 
the p-shell is expected in the wurtzite phase but not in a 



zinc-blende structure. 



B. Single-particle states 

So far we discussed only the role of symmetry on the 
energy spectrum. Turning now to the dipole selection 
rules, one should consider the symmetry properties of the 
wave functions. In the TB description the wave functions 
are expanded in the atomic orbital basis according to 

|</.)=^| S )Ra(R)+^| ft ) R 6 J -(R) , (1) 

R Rj 

where |s)r, denotes the atomic s-orbital centered around 
the site R with (r|s)n = ^> s (r — R), and the correspond- 
ing p-orbitals follow similarly from \pj)u with j — x, y, z. 
The expansion coefficients are given by a(R) and bj(R), 
respectively. 

The action of a symmetry operation T on the wave 
function is defined vi a 52 i 56 

(T^)(r) = ^(T- 1 r). (2) 

In the case of the TB wave function, Eq. ([T|) , this amounts 
to a simultaneous transformation of the orbitals and of 
their centers. Instead of the original orbitals p x and p y , 
it is more convenient to work with the complex combina- 
tions p± = (p x ± ip y )/y2 which are angular-momentum 
eigenstates. 

For a wurtzite QD with C% v symmetry one has to con- 
sider the transformation properties of the orbitals under 
the discrete 27r/3 rotation C 3z around the z-axis and the 
mirror reflection lCi y , defined by (x, y, z) — > (x, — y, z). 
While the atomic s and p z orbitals do not change under 
these transformations, the remaining p± orbitals trans- 
form like 

C 3z \ P+ ) = e-^\p+) , IC 2y \ P+ ) = , 
C 3 ,b_)=e+^b_) , ZC 2 >-> = b+>- 

Using Eq. ([2|) these transformation properties of the local 
orbitals under the group basic rotations are carried over 
to the total wave function, given by Eq. (fTJ). This de- 
fines the action of the group on the TB basis and allows 
the classification of the eigenstates according to the irre- 
ducible representations of the C% v group. States which 
are invariant under the action of the group will be called 
s-states and will be denoted by \<j) s )- This terminology 
and notation are somewhat lax since according to it a 
p z -orbital is also an s-state. Nevertheless we use them 
for the sake of simplicity and in agreement with the lit- 
erature. Similarly, the basis vectors of a two-dimensional 
representation changing under the group operations like 
in Eq. will be called p±-states and will be denoted 
|<£ P± ). 

As an example, it can be shown that an s-state wave 



5 



function must have the form 

10.) = E [ i s >* «( R ) + i**>r ^( R ) ] w 

R 

+ £ [ |p+)r Zr /3(R) + |p_) R 2 R /T (R)] 

R 

where the coefficients a(R), /3 2 (R), and the real part 
of /3(R) are invariant under all rotations of Cs v , while 
the imaginary part of /3(R) is invariant under the proper 
rotation G^ z but changes sign under the action of the 
improper rotation !Ci y . The quantity Zr is defined as 
■^r + '^r or in polar coordinates -J A R + Y^e^ R , where 
and Yr denote the in-plane cartesian coordinates of 
the lattice site R and £r is the polar angle of the point 
R in cylindrical coordinates. It is clear that with a(R), 
/3 2 (R) rotation invariant the first sum in Eq. ^ is s- 
like. In the second sum one has a compensation between 
the p-behavior of the orbitals and that of the Zr coef- 
ficients leading to a combination which is left invariant 
by the rotations. The equivalent of this situation in a 
k ■ p approach would be expressed by the p-character 
of the envelopes multiplying to the p±-type Bloch func- 
tions and compensating the phase factor acquired under 
rotation in order to produce full invariance. In this way, 
the coefficients of the expansion Eq. ((4]) can be seen as 
discretized versions of the envelope function components. 

For the electron and hole wave functions whose modu- 
lus square is depicted in Fig. [Tj the lowest-lying levels 
are indeed invariant under the action of the elements 
of the discrete group and therefore represent s-states. 
As expected, the electron s-state consists mainly of s 
atomic orbitals, i. e. the dominant coefficient is a(R), 
which is C^-invariant as stated above. This state looks 
very similar to what would be expected for an s-state 
in the effective mass approximation: it is symmetric un- 
der rotation, has a single maximum at the center and 
decays to the boundaries. In contrast, the shape of the 
ground state of the holes does not show the same behav- 
ior. It shows the (discrete) rotational symmetry, but has 
a node at the center. This state is expected to consist 
mainly of p atomic orbitals and indeed, an inspection of 
its coefficients shows that the second sum in Eq. Q is 
dominant. Since Fig. [1] displays the sum of the mod- 
ulus square of the expansion coefficients, one obtains 
2\Z^ /3(R)| 2 = 2(X R + y|)|/3(R)| 2 , which explains both 
the rotational symmetry and the vanishing value at the 
origin. In other words the unexpected node at the origin 
stems from the p-character of the coefficients. 

A similar analysis of the coefficients can be carried out 
for the excited states (p-shell), and reveals that these 
states have indeed the \4>p±) symmetry. 

C. Dipole-matrix elements 

The numerical evaluation of the dipole-matrix ele- 
ments shows that the non-zero values of the in-plane 



dipole matrix elements ed^j 1 with e in the x-y-plane are 
much larger than those with e in z-direction. Therefore 
we consider here only the in-plane matrix elements with 
e = l/v / 2( 1,1,0) and denote the corresponding dipole- 
matrix elements by dfj\ Furthermore the matrix ele- 
ments dp* pj , with i,j 6 {x, y} are more than one order 
of magnitude smaller than the other non-vanishing ma- 
trix elements. As a consequence, these contribution can 
safely be neglected in the calculation of optical spectra, 
in which the absolute value of the dipole-matrix elements 
enters even quadratically. As explained below, it turns 
out that these results can be understood on symmetry 
grounds. 

To this end we consider the matrix elements (l z \x ± 
iy\l' z ), from which all the in-plane dipole matrix elements 
can be deduced by linear combinations and are character- 
ized by simple transformation properties. Here we denote 
the states \l z ) by the phase factor e~ llz they pick up un- 
der a rotation around the z-axis of the angle 9, that is, 
l z = for the s-state and l z = ±1 for the p±-states. Us- 
ing the transformation properties of the wave functions 
and those of x ± iy under rotation, one can rewrite the 
matrix elements as 

{l z \x±iy\l' z )=e i ^- l '^ e {l z \x±iy\l' z ). (5) 

In a Cqou group the angle 9 is arbitrary and one would 
obtain the familiar result that the dipole matrix elements 
df 1 v can only be non-zero for l z — l' z = ±1. In particular 

one has df^ = = d p ^ p _ with i, j £ {±1} and only d e s h p . 
and df, h „ can be non-zero. 

In case of the C$ v symmetry the allowed values for the 
angle 9 are only the integer multiples of 27r/3. Then one 
finds the weaker result that the condition for non-zero 
dipole matrix elements is l z — l' z = ±1 modulo 3. This 
means that the discussed symmetry reduction opens ad- 
ditional decay channels. One still has d e Jl = d p h +p+ = 

dp h p _ — but, in addition to dp^ ,,, d^ h p± , now also the 
matrix elements d p f ^_ p _ and d e p h __ p can be non-zero. Nev- 
ertheless, since the last two are vanishing in case of a 
continuous rotation symmetry, it is quite plausible that 
they remain small even in the case of a three-fold axis. 
This behavior is found by the numerical evaluation of the 
dipole matrix elements. In our example of the small QD 
we find \df +p _ | 2 = 4.84 x 10~ 4 (A) 2 , which is negligible 
in comparison with |dfp ± | 2 = 4.73 (A) 2 . 

As we will see in Section |V] these selection rules, and in 
particular the vanishing dipole matrix element d e J£, have 
very important consequences for the optical properties of 
nitride QDs. 



D. Coulomb matrix elements 

It is possible to use the same kind of symmetry ar- 
guments for the calculation of the Coulomb matrix el- 
ements. Again we find that a three-fold symmetry axis 
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TABLE I: All nonzero electron-electron Coulomb matrix ele- 
ments determined numerically from the TB-wave functions. 
The index denotes the s-state \<j> a ) and ± denote the states 
\4>p±}- The matrix elements below the horizontal double 
line would be zero in the case of a Coov symmetry. Explicit 
values are given for of the small the small QD but without 
the inclusion of the internal electrostatic field. 





(i,j,k,l) 




VgU [meV] 




(0,0,0,0) 




93.8459 


(o,+,+,o) 


, (0,-,-,0) , (+,0,0,+) 


, (-,o,o,-) 


81.6389 


(-,-,-,-) , (• 


f ,+,+,+) , (+,-,-,+) , 


(-,+.+,-) 


75.8542 


(0,0,+,-) , 


(0,0,-,+) , (+,-,0,0) , 


(-,+,0,0) 


17.0949 


(o,+,o,+) 


, (+,0,+,0) , (0,-,0,-) 


, (-,o,-,o) 


17.0949 




(+,-,+,-) , (-,+,-,+) 




9.3997 


(+,o,-,-) , 


(+,+,0,-) , (-,0,+,+) 


, (-,-,o,+) 


-0.1972 


(0,+,-,-) , 


(+,+,-,0) , (0,-,+,+) 


, (-,-,+,o) 


-0.1972 



allows for more non-zero elements than a Coo Z -axis. Nev- 
ertheless, these additional matrix elements are rather 
small. If the entire system is rotated by an angle 9 = ^ , 
each single-particle wave function acquires a phase factor, 
while the distance |r — r'| is not affected by this rotation. 
Therefore we find for the interaction matrix elements: 

V ijkl =e^ + ^i-^V ijkl . (6) 

From this one can deduce that Vijki must be zero if 
Af + A| — A| — \\ ^ modulo 3. The z-projection of the 
angular momentum is conserved only modulo 3, and not 
exactly, as it would be in the case for a continuous rota- 
tion axis Cocz^. This gives rises to additional non-zero 
matrix elements V^ki ■ Examples for the electron-electron 
V^ki are given in Tab. U using here the (s,p+,p_) repre- 
sentation. Clearly these additional matrix elements are 
rather small. This is in accordance with our expectation 
about matrix elements that would altogether vanish in a 
system with a higher symmetry. 

Symmetry is also responsible for some Coulomb matrix 
elements having equal values. Such cases are grouped in 
Tab. Q] between horizontal lines. 



IV. MANY BODY PROBLEM 
A. Hamiltonian 

So far we discussed only the single-particle properties. 
However, the investigation of optical properties is an in- 
herent many-particle problem. The Hamiltonian H that 
describes the system of interacting electrons and holes in 



a QD consists of two parts: 

i<j ia 

ijkl 

-J2 V ^%i h l e l' e ^h la . (7) 

The free part Hq contains information about the single- 
particle spectrum e( e ' /l - 1 and describes a system of non- 
interacting charge carriers. Here ei„ (ej^) are annihila- 
tion (creation) operators of electrons with spin a in the 
one-particle states |z) of energy ef. The corresponding 
operators and single-particle energies for holes are hi a 
{h\ a ) and e^, respectively. The electron-electron, hole- 
hole and electron-hole Coulomb interaction between the 
charged carriers is included in Hq ou i . The Coulomb ma- 
trix elements Vij t ki are determined from the TB single 
particle wave functions, as described in Sec. Ill CI 

B. Configuration interaction (CI) 

In a semiconductor QD the finite height of the con- 
finement potential leads to a finite number of localized 
states as well as to a continuum of energetically higher 
delocalized states. If one considers only the energetically 
lowest single-particle states, the eigenvalue problem for 
a given number of electrons and holes has a finite (al- 
beit large) dimension and can be solved without further 
approximations. 

As the Hamiltonian H conserves the total number of 
electrons N e and holes Nh, the Hamiltonian matrix falls 
into subblocks with basis states corresponding to uncor- 
rclatcd many-particle states of the form 

= nc e lr ? Ilc^)" j - 1°> > ( 8 ) 

i 3 

with the occupancy numbers n e / L — 0, 1 and J^. n\ = N e , 
J"V = Nh ■ In order to find the eigenvalues and eigen- 
functions of the interacting problem, the Hamiltonian H, 
Eq. |7|). is expressed in terms of the uncorrelated basis 
states \4>), Eq. ([5]). The resulting Hamiltonian matrix is 
then diagonalized numerically. In this way one finds an 
expansion of the interacting eigenstates for a given num- 
ber of electrons and holes in terms of the uncorrelated 
basis states of the system. These states can be used to 
calculate, for example, the interband emission spectra 
between the interacting Coulomb correlated eigenstates 
of the QD system using Fermi's golden rule: 

J M = fEl(^I^IWI 2 ^-i?/-M . (9) 
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Here \4>i) denotes the correlated initial state with energy 
Ei and \<j>f) and Ef the corresponding quantities of the 
final states. A similar equation holds for the absorption 
spectrum. The Hamiltonian H<\ describes the light mat- 
ter interaction in dipole approximation 



Ha 



-e|E| 



n re„ m 



(10) 



with E being the electric field at the position of the QD, 
e p is the polarization vector, and e is the elementary 
charge. Furthermore the states \n) and \m) denote hole 
and electron single-particle states, respectively. Fermi's 
golden rule, Eq. ((9]), implies that the optical field always 
creates or destroys electron-hole pairs. Hence the initial 
and final states differ by exactly one electron- hole pair. 
For the following discussion it is convenient to express the 
dipole Hamiltonian in terms of the interband polarization 
operator V according to = — elEKT 5 + V'). 

If the analysis is restricted to s- and p-shells and the 
small dipole matrix elements d p h p+ and d^_ p _ are ne- 
glected, one can split V into Pi ow and "Phigh according 
to: 



T^high — ^ ^ {^ps ^p+jCT^s,— & ~t~ ^ps &P-,crh, 



S,—<7 I 1 



V^ —S^ (d eh e h + d eh * e h 



(11) 



This is motivated by the fact that the single-particle en- 
ergy separation for the electrons is larger than for the 
holes. 



V. OPTICAL PROPERTIES 



h(o/ meV 

FIG. 2: Absorption spectrum of the initially empty nitride 
QD. The high-energy side stems from a transition in which 
mainly the s-shell of the hole and the p-shell for the electrons 
is occupied. For the low-energy side the s-shell of the electrons 
and the p-shell of the holes yield the main contribution. 



of the large energy splitting between the s- and the p- 
shell of the electrons, the ps-excitonic transition is well 
separated from the sp transition. Without Coulomb in- 
teraction, these transitions can be found at the sum of 
the single-particle energies e p + and e e s + e p , respec- 
tively. For the absorption spectrum the unusual selection 
rules described in Sec. IIII Cl and the resulting "skew" ex- 
citons lead only to quantitative changes. Even if the op- 
tical spectra could be described in terms of diagonal ex- 
citons, one would still obtain two lines in the absorption 
spectrum^, one from the s-shell and one involving the 
p-shell carriers. However, for the emission spectrum the 
changes are quite important. Most importantily the ex- 
citon and biexciton ground state for small InN/GaN QDs 
remain dark^2S This is the case because the excitonic and 
biexcitonic ground states are dominated by those config- 
urations where all the carriers are in their energetically 
lowest shell together with the fact that the dipole-matrix 
element d e Jl involving these states vanishes. 



To be able to distinguish between effects stemming 
from the symmetry group of the QD and those addi- 
tionally introduced by the internal electric field, we will 
first present the results where the built-in field is artifi- 
cially switched off. In a second step, we will discuss the 
influence of this field on the multi-exciton spectra. 



A. Excitonic and biexcitonic properties 

The considered small QD confines only three bound 
electron states (s- and p-shell) for each spin polarization. 
These states enter the CI approach, together with the 
hole s- and p-shells, which are energetically well sepa- 
rated from the other localized states. The corresponding 
excitonic absorption spectrum shows the two lines de- 
picted in Fig. [21 Here we denote an electron-hole pair, in 
which the electron has mainly a and the hole mainly (3 
character, as an cv/3-exciton with a, (3 £ {s,p}- We found 
in Section IIII CI that in the nitride case the transitions 
do not originate from 'diagonal' ss- and pp-excitons, as 
one may expect, but from sp- and ps-excitons. Because 



B. CI results for multi-exciton emission spectra 

Only for more than two excitons the CI calculation 
provides significant population of the electron and hole 
p-shells. Then emission processes involving the 'skew' 
excitons can take place. As mentioned above the low- 
and high-energy side of the spectrum can be attributed 
to processes where an s-electron or a p-electron is dom- 
inantly involved in the recombination process, respec- 
tively. A schematic representation of the level-structure 
and electron-hole pairs typically involved in high- and 
low-energy transitions is depicted in Fig. 

An inspection of the emission-spectra reveals a blue- 
shift as the number of excitons is increased. This is in 
strong contrast to the results known from the InGaAs 
systemi££L££ but can be explained in terms of the diag- 
onal Hamiltonian, discussed in the next section, and the 
fact that the envelopes for the electrons and holes differ 
strongly. Also it can be seen that, with the exception of 
the bX — > 4X transition, all spectra are rather similar 
if one compares the line structure of the low- and high- 



FIG. 3:(color online) For the example of a three exciton state, 
the light (red) and dark (blue) shaded areas connect the car- 
riers that will lead to an emission at the low- and high-energy 
side of the spectrum, respectively. 

energy side. On the other hand, the oscillator strengths 
of the peaks on the high-energy side are systematically 
weaker than the corresponding ones at the low-energy 
side. These aspects will be addressed in the following in 
detail. 



C. Diagonal approximation 

As discussed in Ref. l57l an approximate description us- 
ing a Hamiltonian that is diagonal in the free states can 
be motivated by inspecting the relative importance of the 
various Coulomb matrix element. The eigenvalues of this 
approximate Hamiltonian are then the non-selfconsistent 
Hartree-Fock energies. This approximation has success- 
fully been used to describe the dominant trends in multi- 
exciton spectra in III-V systems.— £&£2JHL&L 

In a ground state emission spectrum, without Coulomb 
interaction, one would observe one line at the high-energy 
side at Cp+e^ and one at e e s + e p on the low-energy side for 
an initial filling from three to six excitons. With Coulomb 
interaction, one observes instead of single lines clusters 
of peaks as shown in Fig. HJ These clusters blue-shift as 
the number of excitons is increased. The approximate 
position of the clusters can already be explained by con- 
sidering only the diagonal elements of the Hamiltonian 
matrix, see Fig. [5] Note that non-identical envelopes 
lead to additional interaction terms in comparison to 
Refs. 22.57.58. The transition energy AE^^f B from a 
many-particle configuration \A) to a configuration \B) is 
given to a first approximation by££ 

&e,h (12) 

Here e and h denote the single-particle states of the elec- 
trons and holes, respectively, that are depopulated in 
the emission process. The index i involves all single- 
particle states except e and h. The quantity D xx stands 



5X->4X 



4X^3X 




2X^ 1X | 

1- 1X-S-0X 

o I 

1150 1200 1250 1300 1350 1400 1450 1500 

hoi/ meV 

FIG. 4: (color online) Ground state emission spectra for a QD 
with different number of excitons. The high-energy side is 
shown in blue (solid lines), the low-energy side in red (dashed 
lines). For the studied system almost no ground state emis- 
sion is observed for exciton and biexciton. As initial states 
the ground states with 3% = are chosen and the internal 
electric field is switched off. 

for the direct Coulomb matrix elements , while the 

exchange matrix elements X xx are given by V^h . Of 
course, these exchange terms contribute only if the spin 
of e or h agrees with the electron or hole state i. The 
first line of Eq. (|T2]) contains the free particle energies e| 
and of the recombining carriers together with the at- 
tractive Coulomb interaction matrix element — D 6 ^; . All 

en 

terms that stem from the interaction of the electron in 
state e with all the other electrons and holes have been 
grouped in the second line. Similarly, the third line con- 
tains the interaction between the hole labeled with h and 
all the other carriers. Explicitly, one obtains for the high- 
energy transition of the 3X configuration: 

pdiag _ e , h _ n eh 

+ 2D% - Df s - D%% - X% (13) 
+ D% + D h £-2Dt s -X% . 

Similarly, one finds from Eq. (TT2")) for the ground state 
transitions on the high-energy side for more than three 
excitons: 

pdiag _ pdiag _r_ A 7T _ Y~ ee _ Y hh 

I2j 4X—>3X ~ ^3X^2X ^-^Hartree ^ pp A S p ' 

E 5X^iX = E AX 6 -^?,X + A-^Hartroe , (14) 

cidiag _ pdiag i A P 

^6X^5X ~ ^5X^iX ' ^-^Hartree ■ 

Note that the Hartree-shift AEHartrec, defined by 
(D e p e p - Df p ) + (Dp!? - Df s ) , would be zero for identi- 
cal envelopes of electrons and holes. The matrix element 
Xpp is given by V p e ^ . The peaks obtained from 
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FIG. 5: (color online) Comparison of QD ground state emis- 
sion spectra using the Cl-approach (solid lines) and the di- 
agonal Hamiltonian (dashed lines) for various multi-exciton 
transitions. In case of six excitons the two results practically 
coincide. All data are calculated in the absence of the piezo 
electric field. 



the diagonal description provide the approximate posi- 
tion of the clusters calculated by the CI, see Fig. [3 In 
particular the smaller shift of the 4X spectrum relative 
to the 3X spectrum as compared to the shifts involving 
more excitons is well described in terms of the exchange 
matrix elements present in the first line of Eq. (fl4|) . The 
corresponding transition energies of the low-energy side 
can be found by changing e <-> h in the equation above. 

While the central position of the clusters is well repro- 
duced by the diagonal treatment, the different splitting 
within each cluster is not explained. A detailed analysis 
of the states involved in the different emission processes 
shows that a semi-analytic description of the spectrum is 
possible. Such a description will be given in the following 
sections. 



VI. SEMI-ANALYTIC DISCUSSION OF THE 
MULTI-EXCITON SPECTRA 




h(o/ meV 

FIG. 6: Ground state emission spectrum for an initial filling 
of three excitons. As initial states the ground states with 
St ot = are chosen. 



A. 3X — > 2X Emission spectrum 

The ground state emission spectrum with an initial 
filling of three excitons, is dominated by pairs of lines 
on the high- and on the low-energy side, see Fig. [6] We 
will analyze only the low-energy side, as the same line of 
arguments can be carried out for the high-energy side. 

Due to the spin symmetry we can restrict ourselves 
to one of the two possible ground states, for which the 
dominant configuration (according to the CI calculation) 
is shown in Fig. [7] A linear combination of configuration 
\A) and \B), enters the 3X-ground state with the same 
amplitude but opposite sign. Nevertheless, one can re- 
strict the discussion to either \A) or \B), because they 
individually produce the same spectrum and possible in- 
terference terms in the expansion of the coupling matrix 
elements appearing in Fermi's golden rule are zero. 

In Fig. [8] configuration \A) is shown together with the 
2X-configuration (in black) that can be reached via the 
action of T^iow Note that the latter is not an eigenstate 
of the Hamiltonian. Instead it will be mixed with other 
states by Coulomb interaction. The predominant contri- 
bution to this mixture is shown in Fig. [8] in light grey. 
Together these two states either form a spin triplet state 
for the electrons and a spin singlet for the holes (ts) or 
a singlet-singlet (ss) state. These ss and ts states are 
split by the exchange Coulomb matrix element and 
can both be observed in the spectrum. The oscillator 
strengths of the corresponding transitions are equal since 
both final states contain the bright 2AT-exciton state with 
the same probability amplitude. Along the same line one 
finds for the high-energy side an approximate splitting 



In the past, optical multi exciton spectra for nitride 
QDs have only been discussed with selection rules carried 
over from the InGaAs system i 27 i 29 i 31 Additionally, the di- 
agonal treatment, which is often suitable for the descrip- 
tion of major trends in the InGaAs syste m 51 i 58 ' 59 i 60 ' 61 , 
can in the present case reproduce only the overall po- 
sition of the clusters, but is by no means sufficient to 
explain the multiplets. Therefore we will discuss in the 
following the different multi-exciton spectra of Fig. [4] in 
detail. 
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FIG. 7: One of the two configurations dominating the three- 
exciton ground states with S' ot = 0. 
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FIG. 8: One part of the dominant configuration of the 3X- 
ground state characterized by S* = +| and = — \, to- 
gether with the 2X configuration (in black) that is created by 
removing the sp-exciton via the action of Vi ov . Additionally 
the most important dark configuration to which the bright 
2X-state couples via Coulomb interaction is shown in light 
grey. 



of 2Xgy which is, however, considerably smaller. Both 
splittings are in good agreement with the CI result and 
explain the dominant peak structure in Fig. [SJ The ra- 
tio of the peak heights on the low- and the high-energy 
side is given in terms of the dipole-matrix elements by 

r ]eh I r ieh 
u sp I a ps ■ 



B. 4X — > 3X Emission spectrum 

The ground state emission spectrum of AX looks sim- 
ilar to the 3X emission spectrum. However, the ener- 
getic splitting within the cluster is larger and the os- 
cillator strength of the two lines of the cluster differs 
from each other. By considering the involved states, 
we find that the difference can be explained by the 
fact that the final states are now doublet-doublet (dd) 
and quadruplet-doublet (qd) states and no longer simple 
ss and ts states as in the 3X case. The correspond- 
ing quadruplet-doublet splitting, which determines the 
energy difference within the cluster, is larger than the 
singlet-triplet splitting. Evaluating the contribution of 
the three different Sl ot = states, yields an approximate 
ratio of 5 to 4 for the different oscillator strengths within 
the cluster and is in good agreement with the CI result. 



ilar block with configurations le s 3e p 2h s 2h p is found for 
the low-energy transitions. Due to this symmetry one 
would expect the same number of lines on the high- and 
on the low-energy side of the spectrum. This expectation 
is fulfilled for all but the 5X — > AX transitions. 

By combining the first four states |V'i,...,4) m Fig- EH 
one can form ss, st, ts and tt spin states. The last two 
states e) m Fig. [Til allow the formation of an st and ss 
state. While the electrons of the first four states always 
occupy both p + and p- , they occupy only the p- state in 
the last two configurations. Therefore one expects that 
the energy of the ss state formed by the states ipi to ip^ 
differs from the energy of the eigenstate created by com- 
bining -05 and -06- By the same token one expects two 
different energies for the two possible si-states. There- 
fore six different energies can be expected. On the high 
energy side of the spectrum four lines are clearly visible 
and another two may be identified on the left side of the 
cluster. For the low-energy side, however, only two lines 
can be observed. 

To obtain further insight, we construct the Hamil- 
tonian matrix generated from the states {|V'i)}f=i an d 
block-diagonalize it by transforming to spin-eigenstates. 
This way one obtains four subblocks: 

Hts =-§, H ' = \[t I) ' (15) 

ffst =-2tl-H„, fl- tt = -2t-| > (16) 

where we introduced the dimensionless parameters t — 
2y/2X^/ 'XH and t = X^/X™ and measured all ener- 
gies in units of X™ relative to £ dia s + 1 + \ with E dia s 
being the energy of the configurations in the diagonal 
approximation. The six eigenvalues of these blocks are 
E ts = -f , £ ss = ±iVl + i 2 , E st = ~2tT 5%/1 + i 2 , and 
E tt — —2t~ |. From these expression one can read off 
that the st and tt spectrum is shifted by —2t relative to 
the ss and ts spectrum, respectively. 

In order to obtain the corresponding oscillator 
strengths of the transitions, one has to calculate the 
matrix elements \ (AX, i|7\igh|5A^, gs)| 2 , where \4X, i) de- 
notes the i-th eigenstate of the four exciton problem and 



C. 5X — > 4X Emission spectrum 

In contrast to the other multi-exciton spectra, Fig. [5] 
shows a clear asymmetry and additional lines at the high- 
energy side of the spectrum. One of the 5X ground states 
is depicted in Fig.QTJ Taking this as the initial state, the 
main contribution to the final states for the high-energy 
side stem from the 6x6 block generated by the con- 
figurations schematically represented in Fig. [TT] These 
states have the lowest non-interacting energy amongst 
those states that can be reached by a removal of one ps- 
exciton from the configurations shown in Fig. 1101 Their 
electronic configuration is given by 2e s 2e p lh s Zh p . A sim- 
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FIG. 9: Ground state emission spectrum for an initial filling of 
five excitons. As initial states the ground states with St ot = 
are chosen. 
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FIG. 10: Dominant ground state configuration for five exci- 
tons with Sl ot = 0. By flipping all spins a second, degenerate 
ground state is produced. 
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\5X, gs) refers to the ground states schematically shown 
in Fig. [T0]of the five exciton problem. Because the eigen- 
state within the one-dimenstional subspaces of the ts- 
and tt-states is not affected by a varied parameter t, the 
oscillator strength does not depend on the value of t. 
In contrast, the heights of the ss- and st-lines depend 
strongly on t as the amplitude of the different states in 
the linear combination varies with t. Denoting the eigen- 
states of H ss with (ai,/3i)' the oscillator strength £ of 
the ss transition is given by 



E = \dl 



eh\2 



_ h 

V2 + 2 



(17) 



A similar analysis can be performed for the st transi- 
tion. In this case one finds 



eh\2 



V2 



(18) 



for the oscillator strength. As H ss and H s t have the same 
eigenvectors and if v = (a, 0) is an eigenvector, so is 
v' = (— (3, a)*, one finds for the si-transitions the same 
dependency of the oscillator strength as for the ss tran- 
sition. For the analysis of the low-energy cluster, only 
the labels e *-> h have to be changed in all the derived 
equations. 

The dependence of the transition energies and the os- 
cillator strengths on the parameter t are depicted in 
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FIG. 11: (color online) Main contribution to the final states of 
the high-energy side of the 5X — > 4X transition with classi- 
fication Sl ot = and, assuming full angular momentum con- 
servation, i* ot = 0. The arrow between electrons or holes indi- 
cates that additional states can be derived from the displayed 
configuration by flipping the spins of the connected carriers 
simultaneously. This way one can create four different states 
from the first configuration and two from the second one. 



FIG. 12: (color online) Energies of the final states involved in 
the high-energy transition measured in units of Xpp relative 
to £ diag + t + i as a function of the dimensionless parameter 
t = 2V2X^ /X£ for fixed t w 0.47. The labeling refers to the 
spin configuration of the final states and whether the states 
are 'bright' or 'dim'. 



Fig. [12] and Fig. [13J respectively. For above discussed 
QD example, one obtains t w —1.2 for the high-energy 
side and t ~ —2.8 for the low energy side. Therefore we 
expect from Fig. [13] that one will only be able to observe 
in the spectrum one of the ss and one of the st lines 
in addition to the ts and tt line. This explains the four 
clearly visible lines in Fig. [9] of the high energy side. In 
contrast, for the low-energy side only two peaks are visi- 
ble. This can be traced back to the fact that one has in 
this case X^ » X™ or t » 2^2. According to Fig. QJ 
this means that the ss and st as well as ts and tt have 
almost identical transition energies. As a consequence of 
this degeneracy only two distinct lines can be observed 
on the low-energy side of the spectrum. 
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FIG. 13: (color online) Oscillator strengths, shown in green 
(solid line) and blue (dashed line), of the two ss transitions 
as a function of the dimensionless parameter t. A clear redis- 
tribution of oscillator strength from one to the other transi- 
tion is visible. The total oscillator strength, red (dotted line), 
however, remains unchanged. The same behavior is found for 
the st transition. 
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FIG. 14: Ground state emission spectrum for an initial filling 
of six excitons. 



D. 6X — * 5X Emission spectrum 

In the case of two shells for each type of carrier, the 
QD is completely filled with six excitons and there is 
only one ground state possible. By removing one sp- 
exciton or one ps-exciton from this configuration, we find 
that the possible final states are exactly those for the 
OX — > IX transitions, only that the occupied sites in 
the OX — > IX problem are now the unoccupied ones. 
Therefore the 6X — > 5X emission spectrum, see Fig. HH 
is very similar to the OX — > IX absorption spectrum, 
with the main difference being that the lines are shifted 
due to the interaction with the 'background' carriers. 



VII. INFLUENCE OF THE BUILT-IN FIELD 

Nitride QDs grown along the c-axis are characterized 
by the presence of a strong internal electrostatic field 
which has a component stemming from the spontaneous 
polarization and a part generated by strain. When this 
field is included in the calculation of the single-particle 
properties, the electron and hole wave functions are spa- 
tially separated from each other which leads to a reduc- 
tion of the oscillator strength i 28 ' 33 i 34 Furthermore the 
single-particle gap and the Coulomb matrix elements are 
altered. The resulting multi-exciton spectra with and 
without the inclusion of the built-in field are compared 
in Fig. [151 The reduced oscillator strength and the over- 
all red-shift of the spectra due to the QCSE are clearly 
visible. In addition to these results of the modified single- 
particle properties, we find a stronger shift of the lines 
with increasing number of excitons in the presence of the 
internal field. This is another result of the strong sepa- 
ration of the electron and hole wave functions which also 
introduces Hartree-shifts in the spectrum. Indeed, if the 
diagonal approximation is applied as outlined in Section 
IV C[ one obtains again the position of the cluster to a 
good approximation. Another difference can be observed 
in the 5X — > 4X spectrum: The number lines changes 
in the presence of the intrinsic field. Since for all other 
excitonic populations the spectra are only altered quan- 
titatively this is a rather surprising result. But this is re- 
solved by noting that the situation without the intrinsic 



field is rather special as one had X^ w X^. Including 
the built-in field leads to a significant deviation of the 
two matrix elements and henceforth to a clear splitting 
of the previously degenerate lines. 

VIII. MULTI-EXCITON EMISSION FOR A 
LARGER QUANTUM DOT 

In order to give a more representative overview, we 
additionally investigate a larger QD with diameter d = 
5.7nm and height h — 2.3nm . It turns out, that in this 
case the energetic order of the two energetically lowest 
hole levels is reversed in the presence of the internal elec- 
trostatic field. A schematic picture of this situation is 
shown in Fig. [TCJ In the absence of the built-in field, one 
still has the 'usual' order with the s-shell being lower 
in energy than the p-shell. Hence, the spectrum looks 
similar to those previously discussed and is here there- 
fore omitted. However, in the presence of the built-in 
field the two-fold degenerate p-shell constitutes the hole 
ground states and has, according to the symmetry con- 
siderations of Section HTT1 non- vanishing dipole matrix el- 
ements with the electron ground state. This is in agree- 
ment with recent k • p calculations^ and experimental 
results for CdSe QDs 63 grown in the wurtzite phase. As 
an immediate consequence the excitonic and biexcitonic 
ground state is bright. However, the corresponding d^l 
dipole-matrix elements are strongly reduced in compari- 
son with the smaller QD due to the stronger separation of 
the electron and hole wave function in this larger struc- 
ture. 

The resulting multi-exciton spectra are shown in 
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FIG. 15: (color online) Ground state emission spectra for the 
smaller QD (d = 4.5 nm and h = 1.6 nm)with (solid lines) 
and without (dashed lines) the inclusion of the internal field. 
Different number of excitons with Sl ot = are chosen as 
initial states. For the studied system almost no ground state 
emission is observed for the exciton and biexciton. 
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FIG. 16: (color online) Schematic representation of a three 
exciton configuration. The light (red) and dark (blue) shaded 
areas connect the carriers that will lead to a emission at the 
low- and high-energy side of the spectrum, respectively. 

Fig. [T7] In contrast to the intuitive picture in which 
first the states with lowest single-particle energy are oc- 
cupied, a strong population of the hole s-shell for more 
than one exciton is found. The ground states are there- 
fore not given by those states with lowest non-interacting 
energy. This is already confirmed by a calculation that 
contains only the Hartree Coulomb terms and can quali- 
tatively be explained as follows: The attraction between 
the electron and hole being in their respective s-shells is 
stronger than the attraction in the case of s- and p-shell 
carriers. Therefore it can compensate the higher single- 
particle energy of the hole in the s-shell. This leads 
already for the biexciton to an occupation of the hole 
s-shell. However, an additional promotion of the other 
hole is not favored in this case because the increase in 
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FIG. 17: (color online) Emission spectra for a larger quantum 
dot (d = 5.7 nm and h = 2.3 nm) with different number of 
excitons and in the presence of the internal field. For the 
emission of each multi-exciton complex the high-energy side is 
shown in blue (solid lines) , the low-energy side in red (dashed 
lines). Due to the reversed level structure for the holes a 
ground state emission is observed for exciton and biexciton. 
As initial states the ground states with Sl ot = are chosen. 

energy due to the stronger repulsive interaction between 



the holes in their s-shell is higher than the energy reduc- 
tion due to the stronger attraction between s-electrons 
and s-holes. As a consequence, the s- and p-shell for 
the holes are equally populated in the biexciton ground 
state. Therefore the biexcitonic line has approximately 
the same oscillator strength as the excitonic one. 

From three excitons on, the ground states are dom- 
inated by configurations in which both the s-shell for 
electrons and the s-shell for holes are fully populated. As 
a consequence, one obtains qualitatively the same spec- 
tra as in the case of the smaller dot with the 'normal' 
order of the shells. However, one finds that the oscillator 
strengths are strongly reduced in the larger system and 
that the Hartree shifts are even more pronounced, due to 
the strong separation of electron and hole wave functions. 



IX. CONCLUSION 

In this work we have investigated the electronic and op- 
tical properties of lens-shaped InN/GaN QDs. Employ- 
ing a tight-binding model where the weak crystal-field 
splitting and spin-orbit coupling for the system studied 
here are neglected, we found an exactly degenerate single- 
particle p-shell. This degeneracy originates from the C$ v - 
symmetry of the underlying wurtzite lattice. This result 
is in particular intriguing in view of recent discussions in 
the literature about the p-shell degeneracy in zinc-blende 
QDs. Based on the microscopically determined single- 
particle wave functions, the dipole and Coulomb matrix 
elements were evaluated. These matrix elements served 
as input parameters for configuration interaction calcu- 
lations and allowed the determination and further anal- 
ysis of optical properties. Our prediction of dark exciton 
and biexciton ground state for small dots is confirmed 
by symmetry considerations. In contrast to other III-V 
material systems, the emission from nitride-based QDs is 
dominated by 'skew' excitons, so that completely differ- 
ent multi-exciton spectra arise. For larger QDs, we found 
that the strong internal electric built-in field can reverse 
the energetic order of the hole states, which results in 
a bright exciton and biexciton ground state. However, 
the oscillator strength is strongly reduced in these struc- 
tures due to the quantum confined Stark effect. By re- 
stricting the analysis to the energetically lowest shells, 
a semi-analytic description of the optical properties was 
possible, leading to a deeper insight into the origin of the 
various emission lines. 
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